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Abstract 

Nonlinear dynamics ( "chaos theory" ) and quantum mechanics are two of 
the scientific triumphs of the 20th century. The former lies at the heart of the 
modern interdisciplinary approach to science, whereas the latter has revolu- 
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tionized physics. Both chaos theory and quantum mechanics have achieved 
a fairly large level of glamour in the eyes of the general public. The study 
of quantum chaos encompasses the application of dynamical systems theory 
in the quantum regime. In the present article, we give a brief review of the 
origin and fundamentals of both quantum mechanics and nonlinear dynam- 
ics. We recount the birth of dynamical systems theory and contrast chaotic 
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motion with integrable motion. We similarly recall the transition from clas- 
sical to quantum mechanics and discuss the origin of the latter. We then 
consider the interplay between nonlinear dynamics and quantum mechanics 
via a classification and explanation of the three types of quantum chaos. We 
include several recent results in this discussion. 

MSC NOS 37K55, 35Q55 

1 Introduction 

"In the beginning, there was Chaos." These words, though somewhat pre- 
tentious, serve as a microcosm for the excitement that has been engendered 
by nonlinear dynamics, which is often called chaos theory among general au- 
diences. Nonlinear dynamics is at the heart of the modern interdisciplinary 
approach to science. Many people, however, appreciate this subject only in 
a very limited sense. Chaos theory was annointed a glamorous field of study 
when James Gleick's introduction to it appeared in print in 1987.c3 Refer- 
ences to chaos are prevalent in popular culture. It has been mentioned, for 
example, in Jurassic Park, the movie Pi, and an episode of The Simpsons. In 
such references, it is often grossly misapplied, demonstrating that although 
chaos is glamorized in popular culture, it is not really understood by the 
general public. People outside the scientific community are far more aware 
of fractals and the pretty pictures that can be created with them than with 
the analysis behind them and the fact that such behavior serves as a good 
model for systems in several scientific fields, including physics, chemistry, 
biology, economics, and geologyHI 



Quantum mechanics has achieved a similar lofty status in the public eye. 
There have been references to it in countless movies, magazines, and televi- 
sion shows. It has even encroached upon the English language. Like dynami- 
cal systems theory, however, it is grossly misunderstood by public audiences. 
The term quantum leap, for example, refers to a very large change, even 
though the quantum regime encompasses quantities that are so tiny that one 
cannot properly analyze the behavior they describe as part of a continuum. 
Indeed, the quantal regime is one of small jumps rather than large ones. 

One of the goals of studying quantum chaos is to combine the paradigms 
encompassed by nonlinear dynamics and quantum mechanics into one co- 
herent theory describing regimes in which both theories are relevant. Unlike 
the two separate concepts, this notion is not well-developed. Additionally, 
quantum chaos is virtually unknown to the general public, despite the fact 
that this subject seeks to reconcile two objects of the world's fascination. In 
scientific circles, the notion of quantum chaos is well-known, but its facets 
are not understood as well as its two underlying theories. The purpose of 
the present paper is to remedy this situation by providing an introduction 
to quantum chaos as well as a brief survey of some of the prevalent ideas in 
this area of research. Toward this end, we introduce some of the fundamental 
concepts of quantum mechanics and nonlinear dynamics before attempting to 
marry these two fields. We contrast chaotic behavior with integrable motion 
and discuss how dynamical systems theory arose from the study of celes- 
tial mechanics. Similarly, we contrast classical and quantum mechanics and 
then recount the origin of the latter. Finally, we present an introduction to 
quantum chaos that includes a classification of its types, a survey of recent 



results, and an attempt to explain what quantum chaotic behavior actually 
represents. 

2 Order and Chaos 

One of the hallmarks of nonlinear dynamics is the concept of equilib- 
ria, which helps characterize a system's behavior-especially its long-term 
motion. There are numerous types of equilibrium behavior that can occur 
in continuous dynamical systems, but such long-time behavior is restricted 
by the number of degrees-of-freedom (that is, by the dimensionality) of the 
system. In other words, one ignores the transient behavior of a dynamical 
system and only considers the limiting behavior as t — > ±00. In dissipa- 
tive systems, one considers the basins of attraction and repulsion of a given 
dynamical system. The extent of the possible complexity of a dynamical 
system's attracting and repelling sets is determined by the dimension of the 
system.t3'E3'Eil Hamiltonian systems (which are conservative) do not possess 
global attractors or repellors, but their dynamics also becomes more complex 
as their dimensionality increases. 

A one-dimensional system may be described by a single (unforced) or- 
dinary differential equation (ODE) of first order. Its phase space is a line. 
All solutions must either approach a steady state or blow up, because the 
topology of the phase line implies that all equilibrium points separate it into 
two distinct regions. Moreover, any observed blow-up must be monotonic; 
there cannot be any spiraling or other complex behavior. 

Systems consisting of either an autonomous pair of first order ODEs or 
a single (unforced) second order ODE have two dimensions. (A forcing term 



corresponds to a non-autonomity, which increases the dimensionality of the 
system when it is suspended into autonomous form£il) Such systems are 
aptly described by a phase plane. Closed trajectories separate the plane into 
two parts, and if the qualitative behavior is different in those two regions of 
phase space, then the trajectories in question are known as separatrices (see 
Figure [l]). Such separatrices are common in Hamiltonian systems, coming in 
the guise of homclinic and heteroclinic orbits. Limiting behavior may include 
steady states, limit cycles, and blow up (which need not be monotonic). Sys- 
tems that are not Hamiltonian may thus exhibit various flavors of attractors 
and repellors. 

The phase space of systems with n G [3, oo) dimensions is embedded 
in W 1 . In addition to the behavior that can show up in systems described 
in spaces of one or two dimensions, those with at least a third dimension 
may exhibit quasiperiodicity, chaotic ("strange") attractors and repellers in 
addition to other manifestations of chaos such as ergodicity. (Once again, 
attractors and repellors cannot occur in Hamiltonian systems, so one must 
distinguish chaotic behavior in those systems from that in dissipative and 
absorptive ones. Such structures may prove to be relevant to the study of 
dissipative quantum chaosBcll) In many contexts, the concept of dimension 
is related to the number of degrees- of- freedom (dof ) of a system. The number 
of dof of a system is defined as the number of variables required to uniquely 
specify its orientation and position in physical spaced This corresponds 
to the number of directions the system may move in configuration space. 
For example, an unconstrained particle in open space may move in three 
different directions. A (holonomic) system with k degrees-of-freedom has a 



2A;-dimensional phase space. We will not treat non-holonomic systemst3E3 in 
the present work. Such systems have velocity spaces of lower dimensionality 
than their configuration spaces so that a k dof nonholonomic system has a 
phase space of dimension n < 2k. Hamiltonian systems are holonomic and 
conservative. They may behave chaotically as well (as long as they possess 
at least two degrees-of- freedom) , although their brand of chaos is somewhat 
different from that in other types of systems. 




Figure 1: A separatrix that occurs in an integrable configuration of a vi- 
brating quantum billiard in a double-well potential. Trajectories inside the 
separatrix behave qualitatively differently from those outside the separatrix. 

Let us now compare periodic, quasiperiodic, and chaotic behavior. To 
contrast the former two, consider the arcade game Asteroids. When the ship 
that the player controls flies off the screen on one side, it appears at the 
corresponding point on the other side, and the same is true of the top of the 
screen with respect to the bottom. In other words, the game's playing field is 
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a 2-torus, which is properly embedded in 3-dimensional space (which is why 
quasi-periodic behavior is possible). More precisely, consider the following 
vector field on the torus: 

x — 1, 

y = u. (1) 

If uo is rational, the line in phase space eventually reaches its initial point 
so that it is periodic, whereas if u is irrational, the line is quasiperiodic, 
approaching every point on the torus arbitrarily closely. Another contrast 
between periodic and quasiperiodic motion occurs with planetary motion. 
According to Kepler's First Law, each planet's orbit is an ellipse with the 
Sun at one of its foci.E-H A better model, however, is one in which planetary 
motion is described by precessional ellipses. The latter motion is a quasiperi- 
odic analog of the elliptical motion that describes the celestial body at any 
given instant. The quasiperiodicity comes from the fact that the properties 
(eccentricity, angle of inclination, etc.) of the ellipse that describe the in- 
stantaneous motion change gradually over time. Interaction with the other 
planets, in fact, leads to this evolution of the orbital parameters. (Directly 
considering such time evolution is a way to incorporate perturbations due to 
the other planets in the solar system as a small perturbation of the two-body 
problem.) 

In addition to the theoretical distinction between periodicity and quasiperi- 
odicity, there is an issue as to whether one can actually observe this difference. 
( Equivalent ly, can nature tell if a number is rational or irrational in this con- 
text?) Every irrational number can be approximated arbitrarily closely as 
a rational one,Ej a fact that is very important for numerical simulations. If 



one considers a single time series in Fourier space, one cannot in principle 
tell the difference between periodic and quasiperiodic motion, because com- 
puters approximate every irrational number as a rational one. (Thus, the 
computer indicates that the result is periodic-though the period might be 
very long.) However, one can distinguish periodicity from quasiperiodicity 
based on the variation of a parameter if one computes multiple time-series 
plots. For periodic motion, one observes that the ratio of given frequencies 
and higher harmonics remains constant, whereas this ratio varies across dif- 
ferent time series in the quasiperiodic case. On a computer, there is really no 
such thing as quasiperiodicity simply because every number is rational. From 
a practical standpoint, however, one can tell the difference between perioic 
and quasiperiodic behavior as long as the reduced form of the rational num- 
ber contains reasonably small integers (that is, as long as the period is not 
too long). One must nevertheless be careful when doing dynamical computer 
simulations, as it may not always be easy to distinguish a quasiperiodic orbit 
and a periodic orbit with a very long period. 

We are now ready to contrast quasiperiodic and chaotic motion. For 
motion to be chaotic, it must satisfy three properties: boundedness, infinite 
recurrence, and sensitive dependence on initial conditions. The first prop- 
erty simply means that one can find a ball of sufficiently large radius that 
contains the chaotic attractor. The second one implies that if one considers 
an arbitrarily small neighborhood about the initial point of a trajectory, it 
will return to the neighborhood infinitely many times. The last property 
means that two trajectories that emanate arbitrarily closely diverge from 
each other at an exponential rate. (That is, the trajectories are character- 



ized by a positive Lyapunov exponent.) Quasiperiodic motion satisfies the 
first two properties but does not satisfy the third. Two nearby quasiperiodic 
trajectories remain "close" to each other in the sense that they only diverge 
linearly. 



Vibrating Billiiirds 




Figure 2: An example of hard Hamiltonian chaos. 

One can also distinguish several different types of chaos. There is Hamil- 
tonian chaos as well as chaotic behavior in dissipative and absorptive systems. 
Both classes of systems exhibit numerous flavors of chaos. In the former, one 
can observe stochasticity (Figure 0), in which a Poincare map displays a 
bounded set within which there is complete disorder. One can also observe 
so-called "soft chaos" or "local chaos" (Figure |^), in which there is some 
fuzziness near otherwise well-defined curves. If the curves are completely 
well-defined in a Poincare map-that is, they lack fuzziness-they describe 
quasiperiodic behavior, and if one has a collection of dots rather than a 
curve, then the depicted motion is periodic. One can also observe chaotic 
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Figure 3: An example of soft Hamiltonian chaos. 

behavior that falls between these two cases. This quasiperiodic route to chaos 
is described by KAM theorySlll 

Finally, one can extend dynamical systems theory to partial differen- 
tial equations (PDEs). A PDE, which may be treated as an infinite set 
of coupled ODEs, has infinitely many degrees-of-freedom.c3 Such systems 
may exhibit chaotic behavior in both spatial and temporal variables. (The 
chaotic behavior discussed above is temporal chaos.) Possible manifestations 
of such spatio-temporal chaos include a spiral wave route to chaos anologous 
to the period doubling route to temporal chaos described by the Feigenbaum 
sequenced Research concerning such spatio-temporal complexity is quite 
active. 

3 The Origin of Dynamical Systems Theory 

The historical evolution of dynamical systems began with the advent of 
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celestial mechanics.!^' E3 The orbits of the planets were described as the 
"Music of the Spheres," and the solar system was treated as something both 
magical and mystical. Considered strong evidence of Divine creation, many 
scientists set out to explain this extraordinary natural symmetry. Nicholas 
Copernicus (1574-1642), an early celestial mechanician, waited until he was 
on his deathbed to publish his heliocentric theory because he knew that it 
would be considered blasphemy. Other pioneers in celestial mechanics, such 
as Galileo Galilei (1564-1642) and Johannes Kepler (1571-1630), also had to 
be careful with that they published on the subject. (Kepler studied celestial 
bodies in part to explain their divine symmetries, whereas Galileo's work was 
treated quite harshly by the Church.) Sir Isaac Newton (1642-1727) used 
Kepler's three laws of planetary motion and three laws of his own to derive 
the inverse square law of gravitational attraction. Essentially, Newton solved 
the (unperturbed) 2-body problem. 

A natural progression of these results was the attempt to derive a solu- 
tion to the n-body problem. According to the intellectual climate prior to 
the 20th century, the universe was a giant orrery that could be completely 
solved. Nobody had yet succeeded in solving the n-body problem, but surely 
somebody would if given the proper motivation. Scientists such as Pierre 
Simon de Laplace (1749-1827), Joseph Louis Lagrange (1736-1813), Simeon 
Denis Poisson (1781-1840), and Spiru Haretu (1851-1912) attempted to an- 
alyze the stability of the solar system by examining the n-body problem, 
but their results were inconclusive. Moreover, Haretu showed that the meth- 
ods they were using were doomed to indeterminacy. The problem was of 
such a magnitude that the development of new methods was essential to its 
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resolution. 

In Volume 7 (1885-86) of Acta Mathematica was an announcement that 
King Oscar II of Sweden and Norway would award a substantial prize and 
medal to the first person to obtain a global general solution to the n-body 
problem. Jules Henri Poincare (1854-1912) had developed new techniques 
for studying differential equations, and he felt that these would provide a 
good intuitive basis for his attempt at this solution. After more than two 
years of study, the nature of the situation began to take shape. One of 
the problem's secrets was revealed by the special case n — 3. Poincare 
proved that there did not exist uniform first integrals other than the one 
that had already been found. This showed that the 3-body problem could not 
be solved quantitatively by Hamiltonian dynamics (by using first integrals 
to reduce the problem to a solvable one of lower dimension), as even the 
restricted three body problem needed two degrees-of-freedom to describe it 
fully. The n-body problem was thus considerably more difficult than anyone 
had realized. Mathematicians would have to change the way they treated 
systems of this sort. They could no longer rely on quantitative methods 
to study the universe, a fact that countered the prevailing philosophy. The 
presence of chaotic behavior showed that determinism did not imply accurate 
prediction, because even a small perturbation of the initial data of a problem 
could cause an arbitrarily large alteration in the behavior of its solution. 

From Poincare's discovery arose dynamical systems theory. In addition to 
showing that the n-body problem could not be solved analytically, Poincare 
discovered the first manifestation of chaotic behavior in the form of homo- 
clinic tangles. Poincare's work served as the foundation for that of several 
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mathematicians and scientists-including such luminaries as George Birkhoff, 
Stephen Smale, Andrei Kolmogorov, Vladimir Arnold, and Jiirgen Moser- 
and the theory and methods he originated now hold prominent places in 
mathematics, science, and popular culture. 

4 Classical and Quantum Physics 

When introducing quantum mechanics to his students, Richard Feynman 
called it "the description of the behavior of matter in all its details and, in 
particular, of the happenings on an atomic scale. "113 Objects on this scale 
behave like neither particles nor waves. The quantum behavior of all atomic 
objects is the same; there are respects in which they behave like particles 
and other respects in which they behave like waves. The reconciliation of 
this particle-wave duality of matter is at the heart of the transition from 
classical to quantum mechanics. 

In quantum mechanics, quantities such as position, momentum, and en- 
ergy play roles as operators as well as variables (depending on the representa- 
tion in use). In classical mechanics, only the latter role is played. Concomi- 
tant with this additional interpretation is the issue of compatibility of ob- 
servable quantities. Incompatible observables do not commute, and so there 
is an 'uncertainty relation' between them.c3 In such relations, the more pre- 
cisely one knows one quantity, the less precisely one can determine the other. 
The canonical Heisenberg Uncertainty Principle expresses this phenomenon 
between position and momentum, a pair of complementary quantities. From 
a mathematical point of view, this example of quantum-mechanical uncer- 
tainty is a consequence of the Fourier Integral Theorem.0 (One must be 
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careful in more general situations, as this result depends on the fact that 
position and momentum are Fourier transforms of each other.) 

Another contrast between classical and quantum physics is that the for- 
mer has a continuous energy spectrum, whereas that of the latter is dis- 
crete. This is perhaps best illustrated by comparing a classical oscillator 
with a quantum-mechanical one. In both cases, one can find 'normal modes' 
('eigenfunctions', 'wavef unctions') ip n (r,t) and their associated eigenvalues 
('eigenenergies'). One then considers an arbitrary superposition of these 
normal modes: 

oo 

^(r,t) = J2 c nMr,t). (2) 

n=0 

Equation ([|) is traditionally interpreted differently in classical and quantum 
mechanics. In the former, ip(r,t) is a linear combination of all the normal 
modes of an oscillator. One can observe this superposition, for example, when 
conducting experiments with a string or a Slinky. One such experiment is 
to demonstrate that 'travelling waves' are a solution of the canonical wave 
equation by showing that they are one example of a superposition of normal 
modes. In quantum mechanics, one can only observe a single wavefunction 
at one time. Instead of being an expression for the degree of expression of 
a given normal mode, the coefficients c n in the eigenfunction expansion for 
ip{r,t) are instead interpreted as a measure of the likelihood that the nth 
wavefunction ip n (r,t) manifests in a given experiment. The superposition 
ip(r,t), then, represents an expectation of possible wavefunctionns rather 
than a linear combinations of observed ones as in classical mechanics. 

In the above discussion, we were careful to indicate that we were refer- 
ring to the discreteness of energy rather than of the eigenvalues some other 
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quantum-mechanical operator, for which one may yet have a continuum. The 
momentum operator, for example, may admit such a continuum. Addition- 
ally, the distinction in the preceeding paragraph is in many ways cosmetic, as 
classical Sturm-Liouville operators have discrete energy spectra." Neverthe- 
less, this distinction is a convenient one to use from a historical and exposi- 
tory perspective-it is extremely useful for elementary discussions of quantum 
mechanics. In truth, however, the situation is subtler and more complicated, 
as there are important dynamical differences depending on whether a quan- 
tum system's spectrum is discrete, continuous, or contains regions with each 
property. In fact, such differences lie at the heart of the search for genuine 
quantum chaos, because any candidate system that might exhibit such be- 
havior must be (spatially) bounded and unforced with a spectrum that is 
not discrete. It must also be fully quantum and describe a finite number of 
particles.lij 

The probabilistic interpretation of quantum mechanics highlights an im- 
portant difference between classical and quantum physics. In quantum- 
mechanical experiments, it is impossible to predict exactly what can happen 
in a given circumstance in the same sense as in classical experiments, as one 
can only observe a single mode at a time. Instead, one predicts the probabil- 
ities of different events, which can then be measured by repeated experimen- 
tation. However, this distinction between classical and quantum mechanics 
is, in some sense, cosmetic. One can, for example, compute probabilities 
in classical mechanics — the key difference is that all observable quantities 
may be measured simultaneously so that the expressions that correspond to 
the off-diagonal elements in a quantum-mechanical matrix calculation must 
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necessarily vanish. This, moreover, is intimately related to what is perhaps 
the fundamental difference between the classical and quantum theories. In 
classical physics, particles can be labeled by their position and velocity at a 
given time and their trajectories are thereby distinguished. In the quantum 
regime, however, particles do not have definite trajectories, so this distinction 
cannot be made. Indeed, in any experiment, one can switch the labeling of 
two identicle particles without altering its outcome. 

Like dynamical systems theory, quantum mechanics has important philo- 
sophical implications.^ The Heisenberg Uncertainty Principle implies that 
making (highly non-perturbative) observations of a phenomenon affects the 
phenomenon itself. This effect, moreover, cannot be minimized arbitrarily 
by altering experiments. There is a minimum disturbance that one simply 
cannot avoid. In classical physics, an observer is important only in a passive 
sense, whereas in quantum mechanics, the effects of an observation can be 
highly nontrivial. 

5 The Origin of Quantum Mechanics 

Now that we have highlighted several distinctions between quantum and 
classical mechanics, let us review its history. The physics community was 
in turmoil as the 19th century faded into oblivion and the 20th century 
began.c3 Amidst this maelstrom lay the origin of quantum mechanics. There 
were many experimental observations that were inexplicable according to the 
firmly grounded classical theory. As with nonlinear dynamics, the theoretical 
answer to these questions required a new way of thinking. At its very core 
lay subjective probability rather than objective determinism. 
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One of the aforementioned experimental observations was that light ex- 
hibits interference fringes and is therefore a wave phenomenon. However, 
erroneous results are obtained if one attempts to explain the photoelectric 
effect using the postulate that light is a wave. One finds that the energy of 
an emitted electron depends only on the frequency of the incident radiation 
rather than on its intensity as one might expect from the classical optics. 

In 1901, Max Planck (1858-1947) observed blackbody radiation, showed 
that energy due to radiation could only exist in the form of discrete packets, 
and introduced his constant h. In formulating his explanation of radiation 
packets, Planck had to abandon the notion that the second law of thermo- 
dynamics was an absolute law of nature. It was instead a statistical law. 
In 1905, Albert Einstein (1879-1955) discovered the photoelectric effect. Six 
years later, Ernest Rutherford (1871-1937) found that an atom has a positive 
central core surrounded by satellite electrons. Such circulating (and hence 
accelerating) particles radiate energy and so — based on classical theory — one 
would expect the electron to collapse into the nucleus. Why, then, does one 
not observe a burst of ultraviolet radiation emitted as an electron spirals into 
the nucleus? Why, moreover, is the frequency spectrum of light emitted from 
an atom discrete rather than continuous? Niels Bohr's (1885-1962) response, 
published in 1913, was his quantum theory of spectra. 

In 1922, Arthur Compton (1892-1962) discovered that photons scattered 
off electrons. Two years later, Wolfgang Pauli (1900-1958) published his 
famous Exclusion Principle, which states that there are no fermion states in 
which two or more particles share the same quantum numbers. In 1925, Louis 
de Broglie (1892-1987) proposed that the wave-particle duality is a universal 
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characteristic of particles in nature rather than just a phenomenon that is 
observed in light £3 He claimed that the wave nature of matter would become 
evident when the magnitude of Planck's constant h could not be ignored. As 
a consequence, one could observe diffraction patterns from beams of particles 
other than photons. 

Because matter exhibits aspects of both particles and waves, one must 
modify the tenets of classical physics. Bohr stressed the need for reconciling 
this wave-particle duality, introducing his concept of "complementarity" in 
1927. In order to accomplish this, quantum theory must account for the 
discreteness of certain physical properties that entered the realm of physics 
before Bohr's atomic model. Normal modes, for example, are a quantized 
phenomenon from classical mechanics. Classically, one considers a super- 
position of such modes to describe the motion of an oscillator, whereas in 
quantum mechanics one considers an expectation of 'normal modes' (that is, 
wavefunctions) to describe the system of interest. 

The previous year, Erwin Schrodinger (1887-1961) wrote an equation 
describing the wavefunction of a particle based on the laws of quantum me- 
chanics he had discovered. c3 This partial differential equation, 

i^m = H ^r,t), (3) 

was similar to the classical equations that describe other types of waves 
such as those describing sound, light, and (especially) heat. When quantum 
mechanics was first postulated, most scientsists studying it spent the ma- 
jority of their effort attempting to solve the Schrodinger equation. Other 
physicists, including Max Born (1882-1970) and Paul Dirac (1902-1984), 
extended quantum mechanics further by incorporating phenomena such as 



spin. 

In 1927, Werner Heisenberg (1901-1976) discovered his Uncertainty Prin- 
ciple, which implies that corresponding to a smaller error in the measurement 
of a particle's momentum must be a larger one in a simultaneous measure- 
ment of the particle's position (and vice versa). For example, if one performs 
an identical experiment many times in which the position of an electron is 
measured (with a given momentum), then measurement of the position does 
not give an identical result for each experiment. A consequence of Heisen- 
berg's principle is that there are compatible variables that can be simulta- 
neously measured and incompatible ones that cannot be. The same year, 
Clinton Davisson (1881-1958) and Lester Germer conducted experiments on 
the wave properties of electrons and thereby demonstrated electron diffrac- 
tion. Also in 1927, Max Born introduced a probabilistic interpretation of 
wavefunctions ip. He claimed that the quantity H^H 2 was properly treated as 
a probability density rather than just as the intensity of a wave. Born's pos- 
tulate is consistent with the results pertaining to the interference of electrons 
and photons. The presence of particles leads to interference fringes, which 
is generated by the wavefunction. At positions for which ||^>|| 2 is large, the 
probability that the particle is found there is concomitantly large. In the 
presence of many particles, their distribution is described by the probability 
density function \\ip || 2 . 

Additional discoveries have molded quantum mechanics into its modern 
form. In 1928, Dirac discovered a relativistic wave equation and predicted 
the existence of the positron. As quantum mechanics continued to develop, 
physicists realized that there were many phenomena not directly encom- 
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passed in the Schrodinger equation, including electron spin and relativistic 
effects. c3 The filling of these gaps led to relativistic quantum mechanics and 
quantum field theories. 

6 Making Sense Out of Quantum Chaos 

When given a mystery to solve, one seeks to find the simplest model 
possible that properly explains the unknown phenomenon. In science as 
in life, it is beneficial to follow the dictums of Occam's Razor. There will 
always be debates among scientists as to whether models contain extraneous 
elements or fail to offer satisfactory explanations, but the basis of simplicity in 
scientific modelling is one that is almost universally followed. A model must 
not only explain the desired phenomena, but it must be accessible to as many 
people as possible subject to that constraint. The success of a model depends 
not only on its relation to reality but also on its practicality. One can, in 
principle, incorporate every intricate detail when modelling phenomena, but 
this serves little purpose if it obscures what is important. Additionally, one 
also strives to incorporate as many phenomena as possible in such a way as 
to provide a useful abstraction of the component of the universe in which one 
is interested. This dichotomy, in fact, is practically an uncertainty principle 
in and of itself. As one complexifies a model, it (in theory) describes reality 
more accurately, but it is also simultaneously more difficult to understand 
and analyze. 

Dynamicists have yet to establish a consensus as to what types of behavior 
constitute quantum chaos. One of the primary goals of the present paper is to 
offer a classification of the types of quantum chaos as well as an explanation of 
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the behavior and analysis that characterizes each type. In the present paper, 
we attempt to extend to the chaotic regime the well-established differences 
between classical and quantum mechanics. Mathematical objects such as 
wavefunction superpositions are interpreted differently in the two subjects 
even when the behavior is integrable, so we generalize such distinctions to the 
case in which a system behaves chaotically. The differences between classical 
and quantum mechanics in the chaotic regime are a natural extension of those 
in the integrable regime, but the behavioral consequences of these differences 
can often be quite profound. 

Dynamical systems theory brought about a revolution in deterministic 
thinking in science just as quantum mechanics fomented a probabilistic rev- 
olution. The study of quantum chaos is an attempt to marry these two ab- 
stract ways of thinking into a coherent whole in order to describe systems for 
which both quantum mechanics and nonlinear dynamics are relevant. More 
specifically, the concept of quantum chaos is an attempt to extend the notions 
of classical Hamiltonian chaos to the quantum regime.^ In the present paper, 
we review a classification scheme for quantum chaos and offer an exposition 
of the types of behavior described by the term 'quantum chaos.'Q There are 
three types of quantum chaotic behavior: "quantized chaos," "semiquantum 
chaos," and genuine "quantum chaos." Our discussion of quantized chaos 
is influenced by those in Chaos in Atomic Physicsu and Chaos in Classical 
and Quantum Mechanics^ and we mostly follow the work of Porter and Li- 
boffH3E5rEl in our discussion of semiquantum chaos. We draw our discussion 
of true quantum chaos from that in Chaos in Classical in Atomic Physics. 
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6.1 Type I: Quantized Chaos 

Quantized chaos, also known as "quantum chaology," is the most fre- 
quently studied form of quantum chaos. This subject is primarily concerned 
with bounded autonomous systems with discrete spectra. Quantized chaos 
concerns the quantization of classically chaotic systems, usually in the semi- 
classical (h — > 0) or high quantum-number regimes. One looks for signa- 
tures of classically chaotic systems on the quantum level. In quantizing a 
chaotic system, one obtains a configuration that though not chaotic in a rig- 
orous sense nevertheless behaves in a intrinsically different manner than an 
integrable system that has been similarly quantized. Nevertheless, the quan- 
tum dynamics of such systems are still affected in a fundamental manner 
by the fact that their classical counterparts are chaotic. It has been shown, 
for example, that all atoms and molecules except the hydrogen atom (and 
related two-body atomic systems) exhibit chaos when treated classically.BE!l 
The quantum dynamics of these systems are not rigorously chaotic, yet their 
associated waves and energies are strongly influenced by the underlying clas- 
sical chaos. The reason a quantum-mechanical system so obtained is not 
rigorously chaotic follows from the discrete nature of its energy spectrum. 
In classical dynamics, chaotic behavior satisfied boundedness, infinite recur- 
rence, and exponential sensitivity. The properties of boundedness and infinite 
recurrence can be applied to the quantum regime in a well-defined manner. 
One might worry that the discrete energy spectrum alters the topology so 
that one has to be careful about what it means for a "trajectory" (a term 
we use loosely in the quantum regime) to recur in an infinitesimal neighbor- 
hood, but this is not actually a problem, because one can consider things in 
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a different space, such as by taking Fourier transforms between the position 
and momentum bases. The concept of exponential sensitity, however, can- 
not be used directly in the quantum regime, because of the so-called "break 
time" phenomenon. In classical dynamics, exponential sensitivity refers to 
the exponentially fast separation of trajectories that began infinitesimally 
close to each other. If one considers Lyapunov exponents in configuration 
space, one obtains a time series that increases linearly at first but then tapers 
off to be roughly constant as "saturation" occursHI In a bounded system, 
two trajectories can only be separated by a finite distance. In tangent spaces, 
however, this restriction is no longer present. (The tangent space at a point 
in a manifold is the set of all tangent vectors at that point Be3) Time-series 
plots of Lyapunov exponents in such spaces thus display a linear increase 
for all time. In the quantum regime, however, one obtains saturation both 
in configuration space and in the tangent space£3 The saturation in the 
tangent space referred to as the quantum break-time phenomenon. It is the 
reason that the concepts of Lyapunov exponents and exponential sensitivity 
break down in quantum physics. Given that the systems under considera- 
tion satisfy boundedness and infinite recurrence rigorously and that there is 
still exponential sensitivity in some sense, it is not unreasonable to call these 
systems chaotic even though they are not rigorously so. Nevertheless, we 
elect not do so, as it is insightful to distinguish this type of behavior from 
that which is rigorously chaotic. The search for so-called genuine quantum 
chaos is equivalent to the search for a fully quantized system that is chaotic 
in a rigorous sense. Defining chaos as deterministic randomness may prove 
to be very useful in this effort. (With this definition, a system is nonchaotic 
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if the information contained in a system is logarithmically compressed by the 
algorithms use to compute it. That is, one avoids the use of Lyapunov expo- 
nents by comparing the size of information input to the size of information 
output. )0 

As discussed earlier, the notion of Lyapunov exponents does not directly 
carry over to the quantum regime. Even in the semiclassical or high quan- 
tum number regimes of quantum chaology, one still has problems defining this 
concept. Nevertheless, it is desirable to have a notion of stability that can be 
used for these situations. In classical systems, the Lyapnuov exponent is de- 
fined as the rate at which the largest eigenvalue of a trajectory grows. When 
the Lyapunov exponent of a classical trajectory is positive, the associated 
motion exhibits exponential sensitivity. In computing these exponents, one 
must consider how fast neighboring trajectories spread apart, which leads to 
difficulties as discussed above. (Additionally, one no longer has trajectories in 
the traditional sense.) One may thus define an instability exponent x related 
to the eigenvalues of trajectories near periodic orbits. (They are computed 
in the semiclassical regime for so-called "periodic orbit expansions" and the 
Gutzwiller trace formula.) These instability exponents are the needed gen- 
eralization of Lyapunov exponents. We remark that though we use these 
instability exponents in quantum chaology, there are fundamentally classi- 
cal objects that are based on the concept of periodic orbits. In order to use 
classical periodic orbits in quantum mechanics, one must calculate the action 
integral S over one period and also consider its variation 5S. One uses the 
Bohr correspondence priniple to intepret the expression one obtains for 5S. 
Moreover, when a classically chaotic system is quantized, one obtains scars, 
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which describe the narrow linear regions with enhanced intensity that occur 
in an eigenstate's intensity pattern.E3 

Before we indulge ourselves in quantum chaology, let's consider some 
alternate terminology in the generalization of the notion of classical chaos. 
In so doing, we introduce the notion of 'quantum billiard chaos, 'c3 which 
may also be described in terms of the nodal properties of wavef unctions. 
This notion of chaos corresponds to quantized chaology as it manifests in 
quantum billiards. (As we discuss later, this can occur in billiards for which 
the Helmholtz equation is not globally separable.) It is chaotic in the sense of 
boundedness, infinite recurrence, and the instability exponents just defined. 
We earlier generalized the notion of chaos directly to the quantum regime. 
In this alternate formulation, one retains the classical notion of chaos and 
instead generalizes one if its component conditions-the idea of exponential 
sensitivity. Though this procedure is reasonable, we will instead partition 
quantum chaos into three categories as discussed above, because doing so 
allows us to discuss quantum chaos for a wider class of systems. 

The tools used in the study of quantum chaology include random matrix 
theory, level dynamics, and periodic orbit expansions. The former comes 
into play in considering a system's (Hermitian) Hamiltonian matrix H. One 
defines an "uncorrected" probability density p{H). That is, separate blocks 
of H are uncorrelated so that if H is an n x n matrix, then 

n n 
i=\ j=i 

For the special case in which n — 2, equation (Q) becomes 

P {H) = Pu(#ll)pi 2 (#12W#22)- (5) 
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One can show that 

p(H) = C exp {-A tr(H 2 )) , (6) 

where A and C are constants of integration. Related to this is the Wigner 
distribution 

Pw{x) = -xexp y-jx j . (7) 

One expects such Wignerian statistics to hold for the spectra of complicated 
quantum systems with many degrees-of-freedom, in which the associated 
Hamiltonian is similarly complicated. Wigner statistics Pyy are valid only 
if a system has integral spin and is invariant under an anti-unitary trans- 
formation such as time reversal. (An anti-unitary transformation consists 
of the composition of a unitary transformation with complex conjugation.^) 
Among the appropriate complex quantum systems are atoms and atomic 
nuclei. It has been shown that Wigner's prediction is consistent with experi- 
mentally obtained data describing the spacing of energy levels. The hydrogen 
atom in a strong magnetic field is a "simple" system in which this behavior 
is observed. The underlying chaotic features of the system cause certain sta- 
tistical features of the quantum spectrum to obey predictions from random 
matrix theory.Q 

If each of the anti-unitary symmetries is broken, the nearest neighbor 
statistics are expected to be described by 

„ , x 32x 2 / 4x 2 \ 

Pu{x) = ^exp[ , (8) 



7T 2 \ 7T 

where the subscript '[/' stands for 'unitary' since in the present case the 
system's Hamiltonian is invariant under all unitary transformations. If the 
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system has half-integral spin and retains anti-unitary symmetry, then all 
levels of the system are degenerate in the sense of Kramers£2l The nearest 
neighbor distribution of energy level spacings between degenerate pairs is 
given by 

^') = ^exp(-— j (9) 

where the subscript 'S" stands for 'symplectic'. If all the anti-unitary sym- 
metries are broken, we are again in the case Py. 

The three cases P\y, Py, and P$ are characterized by the Dyson param- 
eter /3, which indicates the degree of level repulsion as x — ► 0. In the three 
types of statistics above, (3 takes the respective values 1, 2, and 4 for the prob- 
ability densities Pw, P\ji an d Pg- The matrix elements of the Hamiltonian 
have a Gaussian distribution and for /3 = 1,2,4 are invariant respectively 
under orthogonal, unitary, and symplectic transformations. The random 
matrix ensembles are given the respective names Gaussian orthogonal en- 
semble (GOE), Gaussian unitary ensemble (GUE), and Gaussian symplectic 
ensemble (GSE). In addition to determining the appropriate level statistics, 
the symmetry class of a given Hamiltonian H also influences physical char- 
acteristics such as a system's localization length (which determines how fast 
wavefunctions decay in a given basis). Additionally, the zeroes of the Rie- 
mann (^-function correspond very closely with the GUE,E3 which leads to 
the well-known connection between random matrix theory and this famous 
function. 

The concept of level dynamics is also useful in studying the quantum- 
mechanical analogs of classically chaotic systems. If a quantum system de- 
pends on an external parameter, such as the strength of an externally applied 
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magnetic field, its energy levels depend on that parameter. One may examine 
the changes in these energy levels as that parameter is adjusted. These so- 
called level dynamics provide useful information about a quantum-mechanical 
system's underlying chaotic structure. One studies the eigenenergies E n of a 
Hermitian Hamiltonian 

H (e) = H + eV (10) 

as a function of the perturbation parameter e. (The functions H and V 
are independent of e.) One identifies the eigenenergies E n with fictitious 
particles. To derive the equations of motion, one starts with the eigenvalue 
equation 

H(e)\n(e)} = E n (e)\n(e)} (11) 

and defines the matrix elements 

V nm (e) = (n(e)\V\m(e)). (12) 

Assuming that H is invariant under time reversal guarantees that eigenstates 
and matrix elements are real. The evolution equations of the level dynamics 
are then 

E n (e) = V nn (e), 

■tr 9 \ nrn 

v nn - z 2_^ E _ E , 

V nm = + ^ V nl v lm + , (13) 

which can also be derived from Hamilton's equations. These equations can 
be solved once the initial conditions at e = are known. Note that the struc- 
ture of equation ([H|) is independent of the specific form of the Hamiltonian 
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H . Moreover, equation (|T3| ) shows that the dynamics of energy levels of all 
Hamiltonian systems can be split into two parts, corresponding to equation 
(|T0|). (This separation of the Hamiltonian into an unperturbed Hamiltonian 
Hq plus a small perturbation eV is reminiscent of Melnikov theory.EJ'tiil) The 
characteristics of a given Hamiltonian are used only via the initial conditions. 
Thus, if the initial conditions are unimportant for a sufficiently large e > e, 
then one can use equilibrium statistical mechanics to compute the statistical 
properties of the system's energy levels. The methods of level dynamics can 
be generalized to those of "resonance dynamics. "Q In equation (0), the in- 
dices n and m range over the number of dimensions of the relevant Hilbert 
space Ti. These equations preserve the Liouville volume, the energy 

n n=£m ' 

and the total coupling strength 

Q = Y,\ V ™\ 2 > ( 15 ) 

which gives infinitely many first integrals (constants of motion), since in- 
creasing the coupling strength e generates an orthogonal transformation in 
7i whose invariants (including the traces of various operators) stay the same. 
Moreover, the dynamical system (113) is completely integrable, so that one 
must understand a system with infinitely degrees-of-freedom and infinitely 
many constants of motion. In order to do this, one could take the point of 
view of statistical mechanics. Consider a stationary distribution 

p = l-/«5-7Q (i 6 ) 

Zj 
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where Z is the partition function (a normalization constant). The 'inverse 
temperature' f3 and 'chemical potential' 7 are determined by prescribing 
mean values E = (E) and Q = (Q). One integrates out the variables V nm 
to give a distribution of the eigenvalues E n One obtains a GOE probability 
distribution in the eigenvalues, which is a bit disturbing because this ensem- 
ble now seems to show up in a context sufficiently general that it may not 
be of much value to the characterization of quantum signatures of classical 
chaos.Bj 

A third tool used in quantum chaology is the study of periodic orbit 
expansions. In order to pursue this field, it was essential to develop semi- 
classical methods that worked in the quantum regime. The first procedure 
to do this was the formalism of periodic orbit quantization.^ The central 
result of this theory is the Gutzwiller "trace formula." In the semiclassical 
approximation, only periodic orbits {p} contribute in the evaluation of the 
level density 



p(E) = tr [8{E - H)] . 



(17) 



One finds that the classical approximation p c {E) of the quantum-mechanical 



traceE3 is 



^>4e 



ih ^ 2 sinh(x/2) 



where E represents energy, 



To(E) 



cxp 



<jgi 

\q\ 



S ,7T 

/ I — — I — 

h 2 



(18) 



(19) 



is the primitive period, x(E) is the instability exponent, S(E) = J pdq (in- 
tegrated over the periodic orbit) is the action integral, and I is the number 
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of times the stable manifold is oriented in the local p-direction. The super- 
position of these smooth classical functions yields an approximate spectrum 
for the quantum-mechanical energy levels. We remark that the above for- 
mulation corresponds to a Feynman path integral approach, so that we are 
integrating in the variable q = {qi,q2,Q3) from some initial point q' to some 
terminal point q" . This formulation has been generalized, but by considering 
a Green's function (that is, a propagatorQ'c3), even generalized versions of 
equation ( fL8|) reduce to a trace. There are several other similar formulas 
in the study of quantized chaos, as discussed by Gutzwiller.o In some sit- 
uations, for example, it is appropriate to have a hyperbolic cosine function 
rather than a hyperbolic sine function. One can also study the Riemann-£ 
function to gain insight into the trace formula. 

There are two ways in which the trace formula can be used to deal with 
quantum-mechanical problems. It can be applied "forwards" to calculate 
the level density of a given quantum system based on purely classical input 
represented by the periods, actions, Lyapnuov exponents, and characteristic 
parameters of the associated classical periodic orbits. It can also be applied 
"backwards" if the level density p(E) is given. In this usage, information 
about the periodic orbits is extracted from the level density by a generalized 
Fourier transform based on the trace formula (since the trace formula is in 
the form of an eigenfunction expansion). The forward application of the trace 
formula is considered more difficult, although it has been used successfully 
on occasion. Gutzwiller, for example, applied the trace formula to a system 
consisting of electrons with an asymmetric mass tensor moving in a Coulomb 
potential,^ an important problem in semiconductor physics. The forward 
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transform is considered difficult for three reasons: The number of periodic 
orbits increases exponentially in a chaotic system, these orbits have to be 
computed numerically, and the trace formula has convergence issues that 
have to be circumvented with appropriate summation prescriptions. 

In the field of atomic physics, there are many systems whose classical 
counterparts are chaotic. Perhaps the most famous one is the rotation of a 
diatomic molecule under the influence of externally applied microwavesBc3 
A polar dimer molecule, such as Csl, is located between two plates of a 
capacitor, which is connected to a pulse generator that periodically charges 
and discharges the capacitor's two electrodes. This process creates a time- 
varying, spatially homogeneous electric field, so that the molecule experiences 
a sequence of electric pulses that couple with its dipole moment. This simple 
situation is a deterministic one, as we are assuming that there are no random 
fluctuations. Consequently, given the initial state of the rotating diatomic 
molecule, one may compute its associated wavefunction for all time. 

The above technique captures the essential physical features of a rotating 
diatomic molecule, but it does not truly explain them. In order to understand 
the dynamics of the present example, one approximates it by restricting its 
rotation to a single plane (rather than three-dimensional space) and by ig- 
noring the motion of its center of mass. In this approximation, one treats the 
rotating diatomic molecule as a kicked rotor, an example that has become a 
paradigm of quantum chaology. It has been studied in both a classical and 
quantum setting .BB3 One may treat the classical kicked rotor as a pendu- 
lum (or a rotator) that is agitated at equal time intervals with an impulse 
that varies periodically as a function of the rotor's angular position. The 
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motion of the rotor is then uniform between kicks, each of which changes 
the system's angular momentum. As the time interval between agitations 
becomes smaller, the behavior of the rotor approaches that of an ordinary 
forced pendulum. This classical system behaves chaotically, and the effects of 
this behavior are evident in the behavior of its quantum-mechanical cousin. 
The study of quantum chaology and its applications remains an active 
area of research.^ Theoretical research addresses its analytical structure 
as well as the convergence properties of Gutzwiller's trace formula and its 
derivatives. A major breakthrough occurred in 1986 when Sir Michael Berry 
observed similarities between the trace formula and certain representations 
of Riemann's (^-function. This function has henceforth served as a model for 
studying the analytical properties of semiclassical trace formulas. There are 
more fundamental concerns, however, than the convergence of these formu- 
las. Indeed, it has been argued that there may exist a completely bounded 
chaotic dynamical system whose spectrum (which is both real and discrete) 
is identical to the imaginary parts of the (nontrivial) zeroes the Riemann 
^-function. The existence of this system would accomplish two important 
tasks. As a mathematical problem, it would prove the Riemann conjecture 
that 

C (I + iZ ) = ° (20) 

in the region \Im(z)\ < 1/2 only for z with vanishing real part. As a physical 
system, it would offer important insights regarding the analytical connection 
between classical chaos and quantum energy levels. Research in applications 
of quantized chaos include attempts to use the semiclassical methods in this 
area as a mathematical tool for studying classically chaotic systems. These 
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methods are currently thought to be much more useful for interpreting quan- 
tum spectra and wavefunctions than they are for accurately predicting the 
spectra of classically chaotic systems.Q 

There is active research in other areas of quantum chaology as well. For 
example, there have been several recent discoveries concerning the connection 
between quantized chaos and the three matrix ensembles discussed earlier. 
Additionally, it is not completely understood why the statistics of random 
matrix ensembles are so accurate in describing how classical chaos induces 
universal fluctuations in energy levels. Other topics of current interest in- 
clude diffraction and refraction corrections in semiclassical procedures. Some 
scientists are also studying whether the presence of chaos in a system can 
increase tunnelling. In this context, a wave would tunnel between two is- 
lands in a chaotic sea. It has been surmised that the presence of chaos would 
increase the tunnel splitting of energy levels by several orders of magnitude. 
Lastly, quantum chaology in dimensions greater than two is virtually unex- 
plored. The methods in this subject are only expected to be accurate to order 
h 2 independent of dimension, although this has not been shown rigorously. 
The extension of quantum chaology to higher dimensions should nevertheless 
prove quite fruitful. 

6.2 Type II: Semiquantum Chaos 

Semiquantum chaos concerns systems with both classical and quantum 
subsystems. It can arise, for example, in the form of the dynamic Born- 
Oppenheimer approximation,^ which shows up commonly in the study of 
mesoscopic and chemical physics. This adiabatic approximation arises natu- 
rally in systems that may be expressed as the coupling of slow and fast sub- 
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systems. In the Born-Oppenheimer scheme, chaos may occur in both the clas- 
sical (slow) and quantum (fast) subsystems, although considered separately, 
neither of those regimes is necessarily chaotic. The first step in the Born- 
Oppenheimer approximation is to quantize the fast (electronic) subsystem. 
The second step is to quantize the nuclear (slow) subsystem. If, however, 
the electronic energy levels are too close together, the Born-Oppenheimer 
approximation breaks down, as the electronic and nuclear subsystems are 
nonadiabatically coupled. When this occurs, one treats the nuclear degrees- 
of-freedom as classical variables, thereby obtaining a semiquantal regime in 
which a classical system is coupled nonadiabatically to a quantum one. It is 
this nonadiabatic coupling that produces semiquantum chaos in the resulting 
system, which is sometimes called a semiclassical quantizationM'^B 

One may mathematically abstract numerous systems that exhibit semi- 
quantum chaos as quantum billiards with vibrating boundariesBBE3'&El'Ej 
Such systems are not necessarily expressible precisely as vibrating quantum 
billiards, but such billiards serve as a useful toy model in that that they cap- 
ture many of the features of molecular systems. It is in this abstract context 
that we discuss the notion of semiquantum chaos. 

Quantum billiards describe the motion of a point particle of mass m un- 
dergoing perfectly elastic conditions in a bounded domain of mass M 3> rriQ 
in a potential V. The particle's motion is described by the Schrodinger 
equation with Dirichlet boundary conditions. One defines the "degree-of- 
vibration" (dov) of a billiard as the number of boundary dimensions that 
vary with time. If the boundary is time-independent, the billiard is said to 
have zero dov. The one-dimensional vibrating billiard and the radially vi- 
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brating spherical billiard have a single dov, and the rectanglular billiard with 
time-varying length and width has two dov. The dov of a quantum billiard 
are its classical ( "nuclear" ) degrees-of- freedom. A zero dov quantum billiard 
exhibits only integrable behavior if it is globally separable.^ (A quantum 
billiard is globally separable if the geometry of the billiard is one in which the 
Helmholtz equation is globally separable.) Two simple ways in which global 
separability is violated are when a quantum billiard has a concave boundary 
component and when a billiard is geometrically composite, although it is be- 
lieved that global separability may be violated in other ways (such as with a 
quantum billiard whose boundary is a quartic ellipse), by analogy with known 
non-integrable classical billiards. When global separability is violated in a 
quantum billiard, the observed behavior is "chaotic" in the sense of quan- 
tum chaology. Billiards with concave boundary components, for example, 
share many of the properties of Anosov diffeomorphisms.Q Composite quan- 
tum billiards such as the stadium billiard (whose boundary consists of two 
semi-circles joined by a pair of straight lines) have also been shown to exhibit 
chaotic behavior £3El In globally separable, zero dov quantum billiards, how- 
ever, one expects to observe primarily quasiperiodic behavior, analogous to 
the Asteroids example discussed earlier. (One must be careful with the term 
"quasiperiodic" in the quantum regime just as one is with the notion of chaos. 
Esseentially, one is looking at the quantization of what was a quasiperiodic 
regime in the classical situation.) The easiest example to visualize is that of a 
zero dov rectangular quantum billiard. Imagine that the boundary acts as a 
mirror, so that there are imaginary billiards adjoining the actual one on each 
of its four sides. Modulo translation, one then recovers the same situation, 
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as we described earlier when we discussed motion on a 2-torus. If a wave hits 
a side of the billiard, the reflection in the mirror behaves just as the motion 
in the Asteroids example. Under perfect reflection, the angle of incidence 
equals the angle of reflection, and so perfectly reflected trajectories have the 
same features as trajectories on a 2-torus with respect to periodicity and 
quasiperiodicity. Although a stationary, globally separable quantum billiard 
is necessarily integrable, we remark that the Toda lattice is an example of a 
dynamical system that is integrable but not separableHlEa 

Consider an s dov quantum billiard. The total Hamiltonian of the system 
is given by 

H(a 1 ,---,a s ,P 1 ,---,P s )=K + J2^ + V, (21) 

where a±, ■ ■ ■ ,a s represent the time- varying boundary components and the 
kinetic energy (which corresponds to the quantum-mechanical Hamiltonian 
of the particle confined within the billiard) is given by 

K = -|V (22) 

A two-term superposition of the nth and qth states (that is, a two-term 
Galerkin projectionist) is given by 

il> nq (x, t) = a n A n (t)i) n (x, t) + a q A q (t)i[) q (x, t), (23) 

where the complex amplitudes Aj(t) are time-dependent because the sys- 
tem has time-dependent boundary conditions. Linear equations with such 
nonlinear boundary conditions thus lend themselves to analysis via Galerkin 
methods just like nonlinear partial differential equations such as the Navier- 
Stokes and nonlinear Schrodinger equations. The present problem is of a 
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type known as a free-boundary problem, in which one does not know a priori 
the shape of the domain.Ej It is well-posed by the specification of Dirichlet 
boundary conditions, an initial radius a(0), and initial momentum -P(O), and 
initial amplitudes A/(0). 

For now, we specialize to the case of one dov, in which only a single 
boundary dimension varies in time. Porter and Liboff have showncil that a 
two-state superposition consisting of the nth and gth states exhibits chaotic 
behavior if and only if the quantum numbers corresponding to stationary 
dimensions of the billiard's boundary are the same in both states. This re- 
sult, which manifests in observed coupling behavior in the electronic states of 
polyatomic molecules,E3H is easily extended to any finite-term superposition 
by considering the terms pairwise. The result then states that there must 
exist one pair of states that satisfies the above condition. In the vibrating 
spherical quantum billiard, for example, a pair of states has to have the same 
angular momentum quantum numbers I and m for chaotic behavior to oc- 
cur. (This result is proven using separability and orthogonality conditions 
of the Helmholz differential operator.el) These conditions are not surpris- 
ing, because the rotational symmetry of the system is invariant under radial 
vibrations. 

In an integrable two-term superposition state of a one dov quantum bil- 
liard, the equations of motion are 

i = ¥- < 24) 



*~£^ 
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where A = 2 (ei|Ci| 2 + ei|C 2 | 2 ), and G\ and C 2 are constants such that |Ci| 2 + 
|C 2 1 2 = 1. (The energy parameter A is necessarily positive because e* > and 
the \Ci\ 2 correspond to probabilities.) A special case of this configuration is 
obtained by considering a single eigenstate. 



Equation ( 25Q has been studied numerically in the context of the bi- 
furcations that can occur when one considers different potentials VxB In 
particular, one observes only saddle-center bifurcations (and generalizations 
thereof). Either all the equilibria are centers or-for sufficiently small energies 
A-some of them are centers and others are saddle points (depending on the 
form of the potential). Saddle connections in this system have been studied 
to some extent using continuation methods, although there is room for quite 
a bit more research in this area. 

In the chaotic case, the equations of motion take the form 

_ _uw _ 2/iPz 
X ~ a 2 Ma' [Zb) 

y = —, (27) 

a z 



Ma ' 



(28) 



P 



a = w (29) 



and 



da a A 

where x, y, and z are Bloch variables,@ a represents a displacement, P is 
its conjugate momentum, M is the mass of the billiard, m < M is the 
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mass of the confined particle, /i > is the coupling coefficient between the 
two eigenstates, V = V(a) is the potential in which the billiard resides, 
lj = (e 2 — ei)/h, e± = (e 2 ± ei)/2, and e\ and e 2 (where e 2 > £\) are the 
energies of the two eigenstates.^ The bifurcation structure of this system of 
equations is a generalization of that observed in the integrable case in the 
sense that only generalized saddle-center bifurcations can occur. As before, 
pairs of stable and unstable directions bifurcate to the center manifold as the 
energy of the system is increased. Additionally, there is evidence of saddle 
connections for this chaotic case. 

For the one- dimensional vibrating quantum billiard, Blumel and his co- 
authorsQ'Q considered a two-term superposition state that they computed 
to have a coupling coefficient fi = 3/4. They computed Liapunov expo- 
nents to show expontial divergence (and hence chaos) in this case. However, 
for given parameter values and initial conditions, one can tell whether the 
configuration is chaotic simply by examining the associated Poincare maps. 
This, of course, does not prove rigorously that the configuration is chaotic. 
Nevertheless, it corresponds to the canonical application of KAM theory to 
engineering and the physical sciences, so we consider this a sufficient demon- 
stration of chaotic behavior in the present context. Figures and [] display 
chaotic behavior in the classical variables, and Figure f| shows chaos in the 
quantum-mechanical Bloch variables. Classical Hamiltonian chaos in the po- 
sitions and momenta lead to quantum-mechanical wave chaos in the normal 
modes, whereas chaos in the Bloch variables corresponds to chaos in the 
quantum probabilities. ell 

The single dov quantum billiard that has been analyzed most extensively 
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Figure 4: An example of chaotic behavior in the Bloch variables in a vibrating 
quantum billiard. 

is the vibrating sphere,E3 although the one-dimensional vibrating billiard was 
studied earlierBB One finds that the coupling coefficient /i depends on the 
geometry of the system (as well as what quantum states one considers), but 
that the general behavior of the system is typified by the fact that only a 
single boundary variable is time-dependent. Even though the vibrations of a 
sphere occur in M 3 , only the radial dimension of the billiard actually depends 
on time. Other single dov quantum billiards include the radially vibrating 
cylinder, the longitudinally vibrating billiard, and the rectanglar billiard in 
which either the length or width (but not both) depends on time. 

One can generalize the above notions to quantum billiards with two or 
more dov. Two examples of such quantum billiards are the rectangular bil- 
liard in which both length and width are permitted to vibrate^ and the 
cylindrical billiard with both radial and longitudinal vibrations. The latter 
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may be useful as a model for carbon nanotubes or quantum wires. The the- 
orems cited above can be generalized, although these generalizations require 
further study. 

Before we proceed to the third type of quantum chaos, it is important to 
discuss what is meant when we say that vibrating billiards exhibit quantum 
chaotic behavior when the coupling coefficient \x is nonzero. What exactly 
constitutes such behavior and why should one care about it? We formulated 
the chaotic configuration of single dov quantum billiards as a five- dimensional 
dynamical system: the displacement a and conjugate momentum P are clas- 
sical variables, and the Bloch variables x, y, and z are quantum-mechanical 
variables-as they are obtained from the probabilities \A n \ 2 and \A q \ 2 . An- 
other way to formulate this system is to use action-angle variables, which 
gives one classical degree of freedom coupled to a single quantum-mechanical 
one.E5l In each interpretation, we treat the dynamical equations as a clas- 
sical subsystem coupled to a quantum one. The chaotic behavior in the 
(a, P)-plane represents wave chaos, as the position a is an argument of the 
individual wavefunctions ip n and ip q . Additionally, chaos on the Bloch sphere 
(equivalently, in the quantum-mechanical action and angle) represents chaos 
in the quantum probabilities. Therefore, both manifestations of the observed 
chaotic behavior have interpretations that go beyond classical Hamiltonian 
chaos. Moreover, the observed chaotic amplitudes along with chaotic waves 
give us a chaotic superposition of chaotic normal modes. This is a hallmark 
of semiquantum chaos. 

Vibrating quantum billiards, though an idealized mathematical model, 
bear import as simple manifestations of semiquantum chaotic behavior. They 
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are important for other-more practical-reasons as well. They may be used 
to describe nonadiabatic coupling in polyatomic molecules. Such behavior is 
relevant, for example, in the study of Jahn- Teller systemsHrKj In the present 
context, chaotic configurations of a one dov quantum billiard are analogous 
to a diatomic molecule with two electronic states (of the same symmetry) 
coupled nonadiabatically by the single internuclear vibrational coordinated 
The former gives the quantum degrees-of-freedom, while the latter produces 
the classical ones. (This situation is easily generalized to ones with three or 
more electronic states. Additionally, the radially vibrating spherical quan- 
tum billiard captures features of particle behavior in the nucleus^ and as a 
simplistic representation of the quantum dot nanostructure.B3 The vibrating 
cylindrical billiard may be used as a model for the quantum wire, another 
microdevice component.^ More importantly, it may also prove useful as 
a model of carbon nanotubes. Additionally, other geometries of vibrating 
quantum billiards may have similar applications in mesoscopic physics. That 
is, they may be used as models of various chemical nanostructures, as they 
describe the nanomechanical (electronic-vibronic) coupling that can occur 
in such devices. u Moreover, vibrating quantum billiards generalize Enrico 
Fermi's bouncing ball model of cosmic ray accelerationBEa Finally, when 
some compounds are placed into liquids, one obtains solvated electrons that 
may be described as oscillating billiard systems£§c3 With such a wide ar- 
ray of possible application, vibrating quantum billiards are a very versatile 
model. They provide a simple illustration of semiquantum chaos, they gen- 
eralize a toy model of Enrico Fermi, and they have already been shown to be 
relevant to areas of nuclear, chemical, and mesoscopic physics. 
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There remains much to be studied about semiquantum chaos in vibrating 
billiard systems. For example, one may analyze Galerkin projections with 
more than two terms (in which the Bloch sphere is generalized) and quantum 
billiards with two and three dov. A quantum billiard with three dov (such as a 
rectangular prism billiard with vibrating length, width, and depth) may also 
exhibit Arnold diffusion and cross-resonance diffusion, because the Hamilto- 
nian (classical) subsystem has a number of degrees-of-freedom equal to the 
number of dovs of the billiard. These problems may also be studied using 
action-angle coordinates, in which the number of degrees-of-freedom of the 
system is more readily apparent. Future work also includes the study of other 
geometries, such as the two dov cylindrical quantum billiard and billiards 
with concentric geometries (such as spheres or disks) whose inner and outer 
radi both oscillate. Note, in particular, that this latter example removes the 
assumption of convexity and may also lead to a generalization of the notion 
of dov. Another effect to incorporate is that of rebound from the collisions of 
particles with the billiard's boundary. Additional work to be done includes 
further studies of bifurcations in vibrating quantum billiards as well as an 
analysis of coupled vibrating quantum billiards, which is important because 
quantum dots are often coupled in arrays of various geometries in laboratory 
settings. It may also be fruitful to extend the analysis of vibrating quan- 
tum billiards to a relativistic setting as well as to studying multiple-particle 
vibrating quantum billiards. Finally, the Galerkin method discussed briefly 
in the present paper may also be useful for analyzing nonlinear Schrodinger 
equations, because linear partial differential equations with nonlinear bound- 
ary conditions are similar in several respects to nonlinear partial differential 
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equations. 

Semiquantum chaos may be considered in other settings as well. The fact 
that such systems exhibit exponential sensitivity in their quantum subsys- 
tems (represented by the fact that the Bloch variables behave chaotically) is 
a hallmark of the traditional notion of chaos in a semiquantal settingBa We 
remark that quantum chaos is a type of 'wave chaos.' This generalization 
is nominal in terms of 'type I' quantum chaos, as it make little difference 
whether a system is quantum-mechanical or classical when one is studying 
spectra or wave manifestations of ray dynamics in classical wave systems in 
areas such as acoustics, optics, and electrodynamics. However, when apply- 
ing this generalization to the semiquantal regime, it is important to note 
that for classical wave systems, having a classical boundary is no longer an 
approximation. This shows that 'type II' wave chaos exists in nature. Such 
waves have been studied in classical electrodynamics. 

It would be fruitful to apply the methods that have been used in the study 
of semiquantum chaos to nonlinear Schrodinger equations, which can be used 
to describe optical waves in certain media as well as superfluid hydrodynam- 
ics.E3 Recall, however, that in quantum mechanics, the classical boundaries 
used in the present analysis are in truth an approximation. The nuclear 
degrees-of-freedom (represented by the oscillating components of the billiard 
boundary) may be quantized, resulting in a purely quantum (though higher- 
dimensional) system. The effect of such a quantization on semiquantum 
chaotic systems has not been completely resolved, although every previous 
attempt has produced an example of quantized chaos. In particular, quan- 
tizing the walls of a vibrating quantum billiard leads to a higher-dimensional 
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system that properly falls under the heading of quantized chaos. 
6.3 Type III: True Quantum Chaos 

Bounded, fully quantized systems that exhibit exponential sensitivity and 
infinite recurrence are genuinely quantum chaotic. (Semiquantum chaos de- 
scribes exponential sensitivity in the semiquantal regime; systems in this 
regime consist of classical subsystems coupled with quantum mechanical ones, 
so they are not fully quantized.) Quantum chaology describes quantum sig- 
natures of classically chaotic systems. This regime is fully quantal, but the 
"chaos" observed cannot exhibit exponential sensitivity, as we discussed ear- 
lier. Indeed, the existence of systems that are quantum chaotic in the above, 
stronger sense remains an open question. From our previous discussion, we 
note that the energy spectrum of such a system cannot be fully discrete. 
Otherwise, such a system could not display exponential sensitivity. 

The following example has been proposed as a possibility of such a quan- 
tum chaotic system. Most scientists do not consider it an example of such, 
although it is certainly a very interesting system. Consider a spin 1/2 particle 
passing through a chain of two different magnets (types A and B), sequenced 
according to the following recursion formula: 

M n+1 = M n o M n _i, M = A, Mi = B, (31) 

where the symbol o denotes the operation of appending one chain of magnets 
to another. For example, M 2 = M 1 o M = BA, M z = M 2 oM x = BAB, and 
M 4 = M 3 o M 2 = BABBA. 

In general, the chain of magnets (which, in principle, can be constructed 
in laboratory settings) induces spin precession. The propagator U of a spin 
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1/2 particle satisfies the recursion relation 

U n +i = U n U n -\ (32) 

and can be parametrized by 

where a x , a y , and a z are the Pauli spin matrices.c3 

We follow Bliimel and Reinhardt and consider the special case in which 
the magnetic field of each of the magnets is aligned along the y-axis. It 
follows that the propagator U n represents a rotation by angle /3 n . This leads 
to to following recurrence relation for (3: 

Pn+i=Pn + Pn-i (mod27r). (34) 

For the initial conditions (3q = Pi = 1, we recover the Fibonacci sequence.^ 
Defining b n = f3 n /2n, one obtains the above recursion relation mod 1, which 
can be written as the map Q : w n i— > w n +i — 

1 1 \ 

(mod 1), (35) 

1 / 

where the vector w = (b n , 6 n _i). The map ( p5| ) is very similar to the Anosov 
mapB C : w n i— > w n+ i = 

/ 1 1 \ 

(36) 

whose chaotic properties are well-known. The map Q shares many of the 
properties of the Anosov diffeomorphism. For example, it possesses a stretch- 
ing direction v\ ' = (l,g) with a corresponding eigenvalue e^ = g > 1 that 
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has a positive Lyapunov exponent. The mapping Q thus exhibits exponen- 
tial sensitivity and chaotic behavior just like the Anosov diffeomorphism. 
The quantum dynamics of spin 1/2 particles in the given magnetic chain are 
thus argued by Bliimel and Reinhardt to be truly chaotic. u The sequence of 
rotation angles /3 n (mod 2tt) is consequently also chaotic. 

If the spin 1/2 particles are prepared in a pure spin state polarized in 
the +z direction, the corresponding occupation probability in the +\z) state 
after the nth section is given by cos 2 (/3 n ), so the population in the +\z) state 
must be chaotic as well. Measuring this occupation probability provides an 
experimental test for the occurrence of quantum chaos in the present system. 

Additional issues are involved, which raises doubt as to whether one 
should consider this system a genuinely quantum chaotic one. For a given 
n, the magnetic chain is not chaotic because of the unitarity of quantum 
mechanics, as there is no exponential instability for fixed n. Instead, chaos 
occurs as a function of the discrete variable n. Moreover, one must con- 
sider the length of the apparatus required to observe the quantum chaotic 
behavior described above. It is well-known that the Fibonacci sequence di- 
verges exponentially, so the number of magnets increases exponentially with 
n. Therefore, the action of this chain of magnets is equivalent to free motion 
of a particle on a ring whose position is measured at the end of exponentially 
growing time intervals. This provides an alternate means of understanding 
this example. Because of the exponential increase in the length of the mag- 
net chain, the "physical flight time" it takes for particles to go through the 
actual apparatus also grows exponentially in n. However, if the magnets are 
exponentially close in n, then the "natural flight time" grows linearly, and 
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the magnet chain would consequently also be chaotic with respect to this 
temporal variable. Tranforming to the rest frame of the moving beam par- 
ticles, the quantum- mechanical description of a spin 1/2 particle traversing 
the chain of magnets is equivalent to the quantum description of a stationary 
spin 1/2 particle perturbed by a sequence of external field pulses. 

Note finally that the dynamics of the case with aligned magnetic fields is 
the only one that has been investigated thus far. Additionally, this system of 
Fibonacci magnets is considered to be a genuinely quantum chaotic by only 
a handful of scientists. Other candidate systems have been proposed,@ but 
the existence of chaotic behavior (in the traditional sense) in fully quantized 
systems remains an open question. Nevertheless, there are some clues con- 
cerning where to look. In order to have a chance at representing the Holy 
Grail of quantum chaos, a system must be (spatially) bounded, finite-particle, 
undriven, and fully quantum with a spectrum that is not discrete.^ 

7 Conclusion 

The fields of nonlinear dynamics and quantum mechanics have both 
achieved their share of attention in popular culture. Each has been a part 
of a scientific revolution. For example, the notion of quantum mechanics 
brought an important probabilistic interpretation to science, and the advent 
of dynamical systems theory showed that determinism did not imply solv- 
ability. The study of quantum chaos, which has been increasingly scrutinized 
in recent years, is an effort to marry these two subjects. In the present pa- 
per, we discussed the historical evolution and some principle ideas of both 
subjects. We then divided quantum chaos into three behavioral subclasses 
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and discussed several examples, methods, and results in each of these areas. 
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Figure Captions 

Figure 1: A separatrix that occurs in an integrable configuration of a 
vibrating quantum billiard in a double-well potential. Trajectories inside the 
separatrix behave qualitatively differently from those outside the separatrix. 

Figure 2: An example of hard Hamiltonian chaos. 

Figure 3: An example of soft Hamiltonian chaos. 

Figure 4: An example of chaotic behavior in the Bloch variables in a 
vibrating quantum billiard. 
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